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We show that the GBW saturation model provides an exact solution to the one-dimensional linear 
transport equation. We also show that it is motivated by the BK equation considered in the satu- 
rated regime when the diffusion and the splitting term in the diffusive approximation are balanced 
by the nonlinear term. 
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1. Introduction 

Perturbative Quantum Chromodynamics (pQCD) at high energies can be formulated in coor- 
dinate space in the dipole picture [jjj]. If we in particular focus on Deep Inelastic Scattering the 
scattering process can be described in this picture as interaction of virtual photon which has just 
enough energy to dissociate into a 'color dipole' with the hadronic target carrying most of the to- 
tal energy. The interaction process is described here by the dipole-nucleus scattering amplitude. 
This amplitude can be modeled and we will focus here on the Golec-Biernat Wiisthoff saturation 
model ^ which includes saturation effects. It was motivated by requirements that at the high 
energy limit of QCD the total cross section for hadronic processes should obey unitarity require- 
ments. At present there are much more sophisticated approaches to introduce these requirements 
in a description of scatterings at high energies 01 [5], ^ [7|, §, ^|. However, one can still ask 
the question if there is any dynamics behind the GBW model or to put it differently is there any 
equation to which formula proposed by Golec-Biernat and Wiisthoff is a solution? And what is the 
role of the initial conditions? In this article we report on answer to these questions provided in [|l(]] . 



2. GBW model and a transport equation 



The GBW amplitude following form GBW cross section and related to it by a (x, r) = 1 J d 2 bN(x, r, b) 
reads (here we are interested in the original formulation): 



N(x,r,b) = d(b Q -b) 



(2.1) 



where b is the impact parameter of the collision defined as distance between center of the proton 
with radius bo and center of a dipole scattering on it, r is a transversal size of the dipole, x is 

the Bjorken variable, Rq{x) = ( ^ J is the so called saturation radius and its inverse defines 
saturation scale, Q s (x) = 1 /Rq(x) and xq, X are free parameters. This amplitude saturates for large 



dipoles r 3> 2Rq and exhibits geometrical scaling which has been confirmed by data [11] 



2.1 Transport equation for unintegrated gluon density 



The dipole amplitude (2.1) can be related to the unintegrated gluon density which convoluted 
with the kj dependent off-shell matrix elements allows to calculate observables in the high energy 
limit of QCD. This relation is the following (after assumption that the dipole is much smaller than 
the target): 

,2,, N e 



f(x,k\b) 



4a s n 2 



k 4 V 2 



d 2 r Mx,r,b) 
— expHk-r)— g— 



where r and k are two-dimensional vectors in transversal plane of the collision and r : 
Performing this transformation we obtain the known result [[T2|] : 



f(x,k 2 ,b) 



27r 2 a s 



d(b - b)Rl(x)k 4 exp [-Rl(x)k 2 



(2.2) 
|r|,Jt=|k| 

(2.3) 
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Now motivated by the fact that function f(x, k 2 , b) exhibits a maximum both as a function of x 
for fixed k 2 and as a function of k 2 for fixed x, we differentiate f(x,k ,b) with respect to x and 
f(x,k 2 ,b)/k 2 with respect to k 2 . We obtain: 



lf(x,k 2 ,b)(l-R 2 (x)k 2 " 



d x f(x,k 2 ,b) = ~J^w-WJ (2 . 4) 



^ 2 = (2 " 5) 

Dividing eqn. (2.4) by (2.5) and rearranging the terms and defining J^(x,k ,b) = f(x,k 2 ,b)/k 2 , 
Y = ln^oA> L = \nk 2 /Ql we obtain: 

d Y ^(Y,L,b) + Xd L ,^(Y,L,b) = (2.6) 

which is the first order linear wave equation also known as the transport equation. As it is linear it 
cannot generate saturation dynamically but it can propagate well the initial condition leading to a 
successful phenomenology [Q]. It describes the change (wave) in the particle distribution flowing 
into and out of the phase space volume with velocity A . This wave propagates in one direction. The 
quantity J?(x,k 2 ,b) gains here the interpretation of a number density of gluons with momentum 
fraction x with the transversal momentum k 2 at distance b from the center of the proton. The general 



solution of ( |2.6| ) can be found by the method of characteristics and is given by: 

&(Y,L,b)=JF (L-XY,b) (2.7) 
One can go back from ( |2T6| ) to ( |2.3| ) using following initial condition at x = Xq: 

^{x=x Q ,k 2 ,b) = ^—d(b Q -b)k 2 exp{-k 2 ) (2.8) 

This initial condition has saturation built in, since the gluon density vanishes for small k 2 . Knowing 
the properties of the linear first order partial differential equation we see that the property of sat- 
uration of GBW was a consequence of the wave solution which relates x and k 2 supplemented by 
initial conditions with saturation built in. We also see that the critical line of the GBW saturation 
model visualizing, the dependence of the saturation scale on x, Q s (x) = go (y) ^ s m ^ act f rom 
the mathematical point of view the characteristics of the transport equation. 

2.2 Transport equation for the dipole amplitude in momentum space 

Similar investigations can be repeated for the momentum space representation of the dipole 
amplitude N(x,r,b) which we denote by (j)(x,k 2 ,b). 

♦k* 2 '*) = /S" expHk " r) ^7^ (2 " 9) 

A nonlinear pQCD evolution equation like the Balitksy-Kovchegov (BK) equation written for 
(in large target approximation) takes quite simple form and can be related directly to the statistical 



formulation of the high energy limit of QCD (see [13] and references therein). Applying this 



transformation to (2.1) we proceed with differentiation similarly as before and we obtain: 

d Y <l>(Y,L,b) + Xd L (l>(Y,L,b) =0 (2.10) 
which is, as before, the transport equation. 
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2.3 Relation to pQCD 



It is tempting to investigate the relation between found transport equation and the high energy 
pQCD evolution equations like [Q, [7|, ||]. Let us focus here in particular on the form of the BK 
equation in large cylindrical target approximation for the dipole amplitude in momentum space 
for which the nonlinear term is just a simple local quadratic expression. The BK equation for the 
dipole amplitude in the momentum space reads: 



d Y ^(Y,k 2 ,b) = ax 



J logk 2 



(j>(Y,k 1 ,b)-a^ 1 (Y,k I ,b) 



(2.11) 



where a = and %(y) = 2y/"(l) — y/"(y) ~ V(l ~~ Y) i s tne characteristic function of the BFKL 
kernel which allows for emission of dipoles and therefore drives the rise of the amplitude. The role 
of the nonlinear term is roughly to allow for multiple scatterings of dipoles which contributes with 
negative sign and slows down the rise of the amplitude. This equation provides unitarization of the 



dipole amplitude [ ]14[ ] for fixed impact parameter and admits traveling wave solution in the diffusion 



approximation! 15]. The analytic solution of (2.11) within the diffusion approximation relying 



on expanding the kernel of (2.11) up to second order and mapping it to the Fisher- Kolmogorov 



equation has been obtained by Munier and Peschanski [16]. It reads: 



<j>(Y,k?,b) = e(b -b)< 
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«r(7c) \Q 2 (Y)J \Q 2 (Y) 
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(2.12) 



where y c = 0.373 and Q 2 {Y) is emergent saturation scale given by: 



Q 2 e -«X\r,)Y-^o S Y-^^=^ + m/Y) 



(2.13) 



By inspection we see that ( |2.12 ) does not obey the transport equation. The problem is caused by 
the diffusion term. However, we can consider the asymptotic regime called "front interior" [fiq], 
region where transverse momenta k is close to the saturation scale Q S (Y) and rapidity Y is large 
and where the condition In 2 



k 2 



Qifj) ) /2(XX"(Yc)Y « 1 i s satisfied. In this regime ( 2.12 ) simplifies 
and after taking derivatives as in the previous sections we obtain the following wave equation: 



dyty (Y, L,b)+ X BK d L (j> (Y, L, b) = 



(2.14) 



where Xbk = d log Q 2 (Y)/dY. In the limit where Xbk does not depends on energy JT7| ] we obtain 



A, 



BK 



-ax'iYc 



(2.15) 



3. Conclusions 



In this note we have shown that the GBW saturation model is the exact solution of a one- 



dimensional linear transport equation of the form ( |2.6| ). We conclude that since ( |2.6| ) is a linear 
equation the saturation property has to be provided in the initial condition. We found that for the 
GBW model this equation is universal for the unintegrated gluon density f{x,k 2 ,b) and the dipole 
amplitude in momentum space <p (x, k 2 , b) but the details of the shape of the wave depends on the 
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initial condition which is different for each of them. We also studied the relation of the transport 
equation to the BK equation in the diffusion approximation. We have shown that in the region of 
phase space where diffusion and splitting processes are of the same order as the nonlinear term, the 
GBW model is consistent with the BK equation. 
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